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Abstract 

Starting from a sequence of independent Wright-Fisher diffusion processes on [0, 1], we con- 
struct a class of reversible infinite dimensional diffusion processes on Aoo := {x e [0, I ] N : 
^2i>i x i = 1} w ith GEM distribution as the reversible measure. Log-Sobolev inequalities are 
established for these diffusions, which lead to the exponential convergence to the corresponding 
reversible measures in the entropy. Extensions are made to a class of measure- valued processes 
over an abstract space S. This provides a reasonable alternative to the Fleming- Viot process 
which does not satisfy the log-Sobolev inequality when S is infinite as observed by W. Stannat 
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1 Introduction 

Population genetics is concerned with the distribution and evolution of gene frequencies in a large 
population at a particular locus. The infinitely-many-neutral-alleles model describes the evolution 
of the gene frequencies under generation independent mutation, and resampling. In statistical 
equilibrium the distribution of gene frequencies is the well known Poisson-Dirichlet distribution 
introduced by Kingman [8]. When a sample of size n genes is selected from a Poisson-Dirichlet 
population, the distribution of the corresponding allelic partition is given explicitly by the Ewens 
sampling formula. This provides an important tool in testing neutrality of a population. 
Let 

oo 



Aoo = {x = (xi,x 2 ,...) G [0,1]* : > x k = 1}, 



k=l 
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and 

oo 

V = {x = ( Xl ,x 2 , ...) G [0, 1] N : xi > x 2 > ■ ■ ■ > 0, x k = 1}. 

k=l 

The Poisson-Dirichlet distribution with parameter 9 > (henceforth PD{9)) is a probability mea- 
sure Tig on V. We use P(#) = (Pi(9), P2(0), ...) to denote the V-valued random variable with 
distribution Tig. The component Pk{9) represents the proportion of the kth. most frequent alleles. 
If u is the individual mutation rate and N is the effective population size, then the parameter 
9 = 4Nu is the population mutation rate. A different way of describing the distribution is through 
the following size-biased sampling. Let Uk,k = 1,2,..., be a sequence of independent, identically 
distributed random variables with common distribution Beta(l,9), and set 

(1.1) X\ = V x ,X e n = (1 - t/i) • ■ ■ (1 - U n -i)U n ,n> 2. 

Clearly (X^,X|,...) is in space Aoo. The law of X^X®,... is called the one parameter GEM 
distribution and is denoted by H|j 6m '. The descending order of Xf,X%, ... has distribution Tig. The 
sequence X k , k = 1,2, ... has the same distribution as the size-biased permutation of n#. 

Let k = 1, ... be a sequence of i.i.d. random variables with common diffusive distribution v 
on [0, 1], i.e., v(x) = for every x in [0, 1]. Set 

oo 

(1-2) Qe,v = Y,Pk(0)5e h . 

k=l 

It is known that the law of @g tU is Dirichlet{6 , v) distribution, and is the reversible distribution 
of the Fleming- Viot process with mutation operator (cf. [2]) 

(1-3) Af(x) = e -j\f(y)-f(x)Hdx). 

For < a < 1, 9 > —a, let {Vk ■ k = 1,2, ...} be a sequence of independent random variables 
such that Vk is a Beta(l — a, 9 + ka) random variable for each k. Set 

(1.4) X 8 1 ' a = V 1 ,X d n ' a = (l-V 1 )---(l-V n . 1 )V n ,n>l. 

The law of X^ a , X^ 01 , ... is called the two-parameter GEM distribution and is denoted by TT 3 ^™ . 
The law of the descending order statistic of X® ,a , X®' , ... is called the two-parameter Poisson- 
Dirichlet distribution (henceforth IIq,^) studied thoroughly in Pitman and Yor |12j . The sequence 
X® k ' a , k = 1,2, ... has the same distribution as the size-biased permutation of Tl a g. It is shown in 
Pitman [10] that the two-parameter Poisson-Dirichlet distribution is the most general distribution 
whose size-biased permutation has the same distribution as the GEM representation (jl.4p . A two- 
parameter "Ewens sampling formula" is obtained in [11] . Let @$ av be defined similarly to Og, u 
with X e k being replaced by X k ' a . We call the law of @e,a,v & Dirichlet(9 , a, v) distribution. 
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The Poisson-Dirichlet distribution and its two-parameter generalization have many similar 
structures including the urn construction in [7] and [3], GEM representation, sampling formula, 
etc.. But we have not seen a stochastic dynamic model similar to the infinitely- many-neutral-alleles 
model and the Fleming- Viot process developed for the two-parameter Poisson-Dirichlet distribution 
and Dirichlet(9 , ot,v) distribution. 

As the first result in this paper, we are able to construct a class of reversible infinite dimensional 
diffusion processes, the GEM processes, so that both Hg em and its two-parameter generalization 
n^ e ™ appear as the reversible measures for appropriate parameters. 

In [J3], the log-Sobolev inequality is studied for the Fleming- Viot process with motion given 
by (jl.3p . It turns out that the log-Sobolev inequality holds only when the type space is finite. 
In the second result of this paper, we will first construct a measure-valued process that has the 
Dirichlet(6 , v) distribution as reversible measure. Then we will establish the log-Sobolev inequality 
for the process. 

The rest of the paper is organized as follows. The GEM processes associated with rig 6 " 1 and 
n^ e ™ are introduced in section 2. Section 3 includes the proof of uniqueness and the log-Sobolev 
inequality of the GEM process. Finally in section 4, the measure-valued process is introduced and 
the corresponding log-Sobolev inequality is established. 

2 GEM Processes 

For any i > 1, let aj, bi be two strictly positive numbers. We assume that 

(2.1) inf h > J. 

Let Xi(t) be the unique strong solution of the stochastic differential equation 



(2.2) dXi{t) = (a; - (ai + bi)Xi(t))dt + y/Xi(t)(l - Xi(t))dBi(t), X, t (0) G [0,1], 

where {Bi(t) : i = 1,2,...} are independent one dimensional Brownian motions. It is known that 
the process Xi(t) is reversible with reversible measure iTa-^ = Beta(2ai,2bi). By direct calculation, 
the scale function of X; L {-) is given by 

Siix) = I" ' h! 



^ 7i /2 y 2ai (i-y) 2b! ' 

By ()2.ip . we have lim x ^i Si(x) = +oo for all i. Thus starting from the interior of [0,1], the 
process X{(t) will not hit the boundary 1 with probability one. Let E = [0, 1) N . The process 

X(t) = (X 1 (t),X 2 (t),...) 
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is then a i?- valued Markov process. Consider the map 

$ : E -> Aoo, x = (xi,x 2 ,...) -> (x), <p 2 (x), ..) 

with 

991 (x) = xi,<p n (x) = x n (l - xi) ■ ■■ (1 - s„-i),n > 2. 

Clearly $ is a bijection and the process Y(i) = <E>(X(t)) is thus a Markov process. Let E 1 := 
[0, 1] N be the closure of E, C{E) denote the set of all continuous function on E, and C%(E) be the 
set of functions in C{E) that have second order continuous derivatives, and depend only on a finite 
number of coordinates. The sets C(E) and C^(E) will be the respective restrictions of C(E) and 
C%i(E) on E. Then the generator of process X(t) is given by 




and can be extended to C^(E). The sets B(E) and B(A (X> ) are bounded measurable functions on 
E and A^, respectively. 

Let a = (ai,a 2 , . . . , ),b = (61,62, . . .), and 

00 

Ma,b = JJ 7I"a fc ,6 fc , = a b = ^ ajb o 
fc=l 

Then we have 

Theorem 2.1 XTie processes X(i) and Y(i) are reversible with respective reversible measures fi a ^ 
and E ab . 

Proof: The reversibility of X(t) follows from the reversibility of each Xi(t). Now for any two /, g 
in B(A OD ), the two functions / o a o $ are in B(E). From the reversibility of X(t), we have for 
any t > 0, 

/ /(y)^ y b(y(t))]H a , b (d y ) = / /($(x))E x [ ($(x(t)))] Aia , b (dx) 

JAco JE 

= / a($(x))E x [/($(x( i )))]/x a , b (dx) 

JE 

g(y)E y [f(y(t))]~ aih (dy). 
Hence Y(t) is reversible with reversible measure H ajb . □ 

Remark. The one parameter GEM distribution Hg em corresponds to a\ = \,bi = |, and the two 
parameter GEM distribution n^ e ™ corresponds to a« = , 6j = . 
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3 Uniqueness and Poincare/Log-Sobolev Inequalities 



Let 

oo 



i=l 

pN 



be the closure of space Aoo in IR under the topology induced by cylindrically continuous functions. 
The probability H a b can be extended to the space Aoo. For simplicity, the same notation is used 
to denote this extended probability measure. 
Now, for x <G Aoo such that 

n 

^ Xj < 1, for all finite n, 



i=l 



let 



oo n2 00 o 

l,J=l J 1=1 

where 



„ m ._ ^ „ (^(1 - Eti 1 x i) - ^XM 1 - Ef=i 1 x i) - x k) 

fc =i x k(l-22i=ixi) 

, / x _ ^ ~ Ef^i 1 gi) - ^fc)(Qfc(l ~ Ef^i 1 x i) - («fc + frfc^fc) 

fc =i W-X,j=i^) 

Here and in what follows, we set X)i=i = and Ili=i = 1 by conventions. By treating jj as one, 
the definition of ££ (x) can be extended to all points in Aoo. Through direct calculation one can 
see that J2? is the generator of the GEM process. 
It follows from direct calculation that 



(3.1) ^2 \aij(x)\ <3, \bi(x)\ < ^~2(b k Xk + a k ), x£A w . 

i,j=l k=l 



Indeed, since 1 — Yl\=\ x i — x i an d Xa<i<j<oo x i x j < \i we obtain 
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^2 \ciij(x)\ = ^2au(x) +2 ^ \aij(x) 

i,j=l 8=1 l<i<j'<oo 



co < 9 ; i— 1 



2 f 1 ~ 2_^=1 x l | ^fe 



Eft 
Z=l^ 



i=i Xl fc=i 1 ~ 

i-l 

+2 £ ^i(i+E r _ v , ,.. 

l<i<j<co fc=l Z^=l^« 

co i i— 1 



i=l 1=1 k=l 

oo co 



1=1 J =2 + 1 

< 1 + 2 = 3. 



Thus, the first inequality in (|3.1j) holds. Similarly, the second inequality also holds. 
Let 

0/0*0 %(a0 



'CO 1 

r(/.#) = E^J fey 



Then r(/,/) G C 6 (A M ) for any / G C£ (£«,). 

For each a > 0, 6 > 0, let a aj b be the largest constant such that for / G C£ ([0, 1]) the log-Sobolev 
inequality 



(3.2) vr a , fe (/ 2 log f 2 ) < — f x(l - x)/ (x) 2 ^^) + vr a , fe (/ 2 ) log ^ a , fe (/ 2 ) 

holds. According to [131 Lemma 2.7], we have a a ^ > Moreover, it is easy to see that for 
a,b > the operator 

Kl-r)^ + (a-(o + 6)r)^ 
on [0, 1] has a spectral gap a + 6 with eigenfunction /i(r) := a — (a + 6)r. So, the Poincare inequality 



I r 1 

(3.3) TT a , b (f) < — — / x(l - x)/'(x) 2 vr aif) (d 2; ) + ^, fe (/) 2 

holds. 
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Let C^([0, 1] N ) denote the set of all bounded, C°° cylindrical functions on [0, 1] N , and 

^ = {f\ Aoc :feC^([0,lf)}. 
Then we have the following theorem. 
Theorem 3.1 For any f,g G ^C£°, we have 

(3-4) *(f,g) := Sa.bOU <?)) = -E^f&g). 

Consequently, <^C£°) is closable in L 2 (A 00 ;H aj b) and its closure is a conservative regular 
Dirichlet form, which satisfies the Poincare inequality 

(3.5) ~ a , b (/ 2 ) < . 1 ,,/ (/, /), / € ^(<T), H a>b (/) = 0. 

If moreover inf{aj A hi : i > 1} > 0, the log-Sobolev inequality 

(3-6) H a , b (/ 2 log/ 2 ) < ^-S{f,f), f G ^(<T), ~ a , b (/ 2 ) = 1 

Pa,b 

/joWs for some /3 a>b > infj^^ : i > 1} > 0. 

Proof: For any f,g£ ^C£°, there exists n > 1 such that 

(3.7) /(x) = f(xi, ■■■ ,x n ), g(x) = g(xi, ■■■ , x n ), x = (x 1 , ■ ■ • , x n , ■ ■ • ) G [0, 1] N . 

Let 

</? (n) (x) = (pi(x),... )¥ > n (x)), 
which maps [0, l] n on to A n := {x G [0, If : £?=i X, < 1}. Define 

n a n 5 

i=l ' i=l 

and 

n 

I I _ i-jti _n _ ,„(n) — ^ 

i=l 

Then, regarding {H n := 7r" b o <p( n ) 1 : n > 1} as probability measures on 



Aoo, by letting ~ n := E n (dxi ■ ■ ■ dx n ) x S (dx n+1 , ■ • • ), it converges weakly to H ab . S inc6 L n is 
symmetric w.r.t. 7r™ b we have 



(3.8) 



/ 0)1]n E-(l--)(^/°^)(^o^)) d<b 
= - / 5 o^)L n /o^)d< b . 



Noting that 



i-i 



(p i (x)=x i T[(l-x l ), Xi = - — J^ff i > 1, 

1=1 !-Ei=i^( x ) 



we have 



Therefore, 



(3.9) 



d/o<pW(x) = (Sjj -Xj)ipj(x) df q (n),, 

dxi 4^ Xi(l - Xi) difj 

J>1 x ' J 



= / r(/, 5 )o^) d < b = / T(f,g)dZ n . 

J [0,l] n J An 



/ [0,l] r 

By ()3.ip and (|3.7p . we have T(f,g) € C^Aoo) so that the weak convergence of H n to H ab implies 

(3.10) Jim/" T(f,g)dZ n = f r(/, 5 )d~ a , b . 

Similarly, by straightforward calculations we find 

L n /o^")(x) = (^/)o^)(x). 
Moreover, (|3. 1 j) and ()3.7|) imply that gJff € C^Aqo). Thus we arrive at 

hm / ff o^)L n /o^)d< b = / 5^/dS a , b . 

Therefore, (|3.4p follows by combining this with (|3.9p and (|3.10p . This implies the closability of 
(tp 5 , J^"C£°), while the regularity of its closure follows from the compactness of under the usual 
metric 
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p(x,y) := ^2 2 l \ x i-ViY 
1=1 

Indeed, it is trivial that 9(g) n C ([0, 1] N ) D ^C fe °° which is dense in 0(<f) under g^ 2 given by 

g(f,f) = £(fj) + 11/111- 

Moreover, for any i 7 € C(A (X) ) = Co (Acq), by its uniform continuity due to the compactness of the 
space, 

AooX i-> F n (x) := F(xi, ■ ■ ■ , x n , 0, 0, • ■ ■ ), n > 1 

is a sequence of continuous cylindric functions converging uniformly to F. Since a cylindric con- 
tinuous function can be uniformly approximated by functions in J^"C£° under the uniform norm, 
it follows that ^C£° is dense in C^Aqo) under the uniform norm. That is, the Dirichlet form 
(g,9(g)) is regular. 

Next, the desired Poincare and log-Sobolev inequalities can be deduced from (|3,3p and (|3,2p 
respectively. For simplicity, we only prove the latter. By the additivity property of the log-Sobolev 
inequality (cf. [5]), 

holds for all h € C#([0, l] n ), where 

/5a,b = infteoi,* : » = 1, • . . , n}, / (n) (x) = . . . , x n , 0, . . .). 
Combining this with (|3.9p . for any f & 9, the domain of .5f, we have 

s n (/ (n)2 log/ („)2 ) < 1 f F (n) {f j )dE n + E n {f (n)* )logE n ( f (nfy 
Pa,b 

Therefore, as explained above, (|3.6h for f & 9 follows immediately by letting n — > oo. Hence, the 

1/2 — 

proof is completed since 9(g) is the closure of 9 under \ 1 . □ 

We remark that since (g, 9(g)) is regular, according to OS], (L, 9) generates a Hunt process 
whose semigroup Pt is unique in L 2 (H a b)- Thus the GEM process constructed in section 2 is the 
unique Feller process generated by Jz? . Moreover, it is well-known that the log-Sobolev inequality 
(|3.6p implies that Pt converges to H a t, exponentially fast in entropy; more precisely (see e.g. [TJ 
Proposition 2.1]), 

2a,b(ft/ log Pt/) < e" 4 ^ b * ~ a , b (/ log /), / > 0, S a , b (/) = 1. 
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Moreover, due to Gross [3], the log-Sobolev inequality is also equivalent to the hypercontractivity 
of Pf 

Thus, according to Theorem 13. 1\ we have constructed a diffusion process which converges to its 
reversible distribution H a b in entropy exponentially fast. 

4 Measure- Valued Process 

It was shown in Stannat [13] that the log-Sobolev inequality fails to hold for the Fleming- Viot 
process with parent independent mutation when there are infinite number of types. In this section, 
we will construct a class of measure-valued processes for which the log-Sobolev inequality holds 
even when the number of types is infinity. 

Let us first consider a measure- valued processes on a Polish space S induced by the above con- 
structed process and a proper Markov process on S' N . More precisely, let X t := (Xi(t), • • • , X n (t), • • • ) 
be the Markov process on Aqo associated to (<^, ^(<^), and £t := (£i(t), • ■ ■ , £ n (t), • • • ) be a Markov 
process on independent of X t . We consider the measure- valued process 

oo 
i=l 

where Xj can be viewed as the proportion of the i-th family in the population, and its type or 
label. Then the above process describes the evolution of all (countably many) families on the space 
S. Let be the set of all probability measures on S. Then the state space of this process is 

:= {7 6 ^1 '■ supp7 contains at most countably many points}, 

which is dense in ^\ under the weak topology. 

Due to Theorem 13. 1\ if £t converges to its unique invariant probability measure v on S N , then 
rjt converges to II := (H a ,b x u ) V 7 " 1 f° r 

00 

ip : x 5 N -» J( Q ; i/)(x, f ) := . 

i=i 

Unfortunately the process r]t is in general non-Markovian. So we like to modify the construction 
by using Dirichlet forms. 

Let v be a probability measure on and (<£gN, &{Sgn)) a conservative symmetric Dirichlet 
form on L 2 (v). We then construct the corresponding quadratic form on L 2 (^q;U) as follows: 

£* (F,G) := f ^,G^(dO+ f ^M(F x ,G x K, b (dx) 

F, G G := {H G L 2 (n) : H x := H o 4>(x, •) € @{£ s n) for S a>b -a.s. x, 

:= Hoip(-,£) G 9{S) for is-a.s. £,such that £* {H, H) < 00}. 
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Since II has full mass on *^o, to make the state space complete one may also consider the above 
defined form a symmetric form on L 2 («^xj H)(= L 2 (^o;TL)). 

Theorem 4.1 Assume there exists a > such that 

K/ 2 log/ 2 ) < -<%"(/,/) + K/ 2 )logK/ 2 ), / G 
a 

holds, then 

(4.1) n(F 2 i og F 2 ) < 1 c# (i ? ,i ? )+n(F 2 )io g n(F 2 ), f g 

a A /3 ajb 

If moreover &{S^g Q ) C L 2 (^#i;II) is dense, then is a conservative Dirichlet form 

on L 2 (^#ojn) s ° that the associated Markov semigroup Pt satisfies 

(4.2) U(P t F log P t F) < IL(F log F)e" (/3a ' bAa) *, i > 0, F > 0, 11(F) = 1, 
and (<C# , £F(<C^ )) is regular provided so is {S s n^{Sgs)) and S is compact. 
Proof: Let 

= {F G L 2 (~ a , b x v) : F(x, •) G ^(<%n) for S a , b -a.s. x, 

F(-,£) G 0(<£) for z/-a.s. £,such that <#(F,F) < oo}, 

where 

i(F,G):= [ ^N(F(x,.),G(x,-))3a,b(dx)+ / ^(F(-, £), £))"(<!£)■ 

Then (<^, is a symmetric Dirichlet form on L 2 (A co x »S N ;H ajb x i/) and (see e.g. [5J Theorem 

2.3]) 

(4.3) (» a , b xi/)(F 2 logF 2 )< 1 (H a , bXi /)(F 2 ), F G ^(<#), (S a , b x z,)(F 2 ) = 1. 

Let Ft be the Markov semigroup associated to (<$", Then (|4.2p follows from the fact that 

rjt = if>(X(t),(;(t)) and implies (cf. pQ Proposition 2.1]) 

(H a , b x j/) (P t G log P*G) < (E a , b x i/)(GlogG)e- 4 ^^Aa) i 

for all t > and nonnegative function G with (H ab x v)(G) = 1. Since F G ^(<C# ) if an d only if 
Fo^ G and 
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g J V {F 1 F)=g(Fa^Foip), 

(|4.ip follows from (|4.3p . By the same reason and noting that ($,!3(E)) is a Dirichlet form, we 
conclude that {S^ x , £F(<C^ )) is a Dirichlet form provided it is densely defined on L 2 (^#i;n). Finally, 
if S is compact then so is ^#1 (under the weak topology). Thus, as explained in the proof of Theorem 
3-H for regular {S s n,2(S s n)) the set 



{/((-, 9i),--- ,(;9n)) : »>l,/e«ji GC(5),1 <i< n }cC„(l )nf(^) 

is dense both in Co(-#i))(= C(.-#i)) under the uniform norm and in Qi{£j(^) under the Sobolev 
norm. □ 

Remark. Obviously, we have a similar assertion for the Poincare inequality: if there exists A > 
such that 

Hf 2 )<\^(f,f) + Hf) 2 , f e °A<$s«) 



holds, then 



n(F 2 )< . A . f -—MF,F) + u(Fy, Fef(^ ) 



To see that the above theorem applies to a class of measure-valued processes on S, we present 
below a concrete condition on <#^n such that assertions in Theorem 14.11 apply. In particular, it is 
the case if S'gn is the Dirichlet form of a particle system without interactions. 

Proposition 4.2 Letvi be the i-th marginal distribution of v and for a function g on S letg^(£) := 
g(S,i),i > 1- Assume that 

^0 ■= \g G C (5) : g {i) G ^(^n), sup# 5N ( 5 », 5 «) < oo) 
is dense in Cq(S). Then (<C# , 3>{$jt )) is a symmetric Dirichlet form. 

Proof: Under the assumption and the fact that CS(Aoo) is dense in L 2 (^q;IV), the set 

& := {/((-, 91 ), ■■■ , (;g n )) : n > 1, / G C£(R»), 5i G 1 < i < n} 

is dense in L 2 (^q;IV). Therefore, by Theorem 14.11 it suffices to show that 5? C £F((f# ); that is, 
for F := f((;gi), • • • , (;g n )) G Y, one has F o r/> G 0(<f). Let 

F m (x) = i^^a^i (&),••• ,^£iSn(£i)J, x G Aoo, m > 1. 
i=i i=i 
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Since for fixed £ G S 1 , 

n 

fi^Fo^^x)^^), i>l 
fc=i 

is uniformly bounded, one has F m G £F(<f ) and (|3.1|) yields 

&(F m , Fm) < C 

for some constant C > and all m > 1 and £ € S" N . Thus, F o G 3>(£) for each £ G S N and 

(4.4) sup o </>(•, £),Fo </>(•,£)) <C. 

{ 

On the other hand, since gt G 1 < A; < n, noting that for any x G 

n oo 

|F o V.(x, i) — F o ^(x, £')| 2 < (E ll^/Hoo) 2 E " »fc(^i)| 2 . 

fc=l 8=1 

we conclude in the spirit of [9, Proposition 1-4.10] that F o ?p(x, •) G ^(<^n) and 

^(Fo^(x,),Fo^(x,.))<C' 

for some C" > independent of x. Combining this with (14.411 we obtain F o ^ G £F(<f ). □ 
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